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ABSTRACT 

The Virgo and A1795 clusters of galaxies were re-observed by EUVE with in 
situ background measurements by pointing at small offsets. Earlier, a similar re- 
observational strategy applied to the cluster A2199 revealed that the background 
radial profile was consistent with a flat distribution, and therefore the original 
method of extracting cluster EUV signals by the subtraction of an aymptotically 
determined background was valid. It is shown here that the same conclusions 
hold for the current sample. A model of the background was obtained from its 
known properties and the in situ measurements, and the subtracted cluster fluxes 
remain in agreement with those reported in our discovery papers. They are also 
consistent with results from the most conservative procedure of direct point-to- 
point subtraction of the in situ background and proper error propagation, which 
still preserves the existence of the EUV excess and its rising radial trend. We 
present evidence which argues against the soft excess as due to peculiarities in 
the line-of-sight Galactic absorption. The data appear to favor a thermal origin 
of the emission. 

1. Introduction 

Clusters of galaxies are sources which emit EUV and soft X-ray radiation (~ 0.1- 0.4 
keV) substantially in excess of that expected from the hot intra-cluster medium (ICM) 
at X-ray temperatures (a few keV). The Extreme Ultraviolet Explorer mission (Bowyer 
and Malina 1991) was pivotal to the discovery of this phenomenon, as it opened a unique 
window to the EUV sky with the Deep Survey (DS) Lex/B filter (~ 65-190 eV, hereafter 
the EUV band). Cluster soft-excess (CSE) emission was reported for the Virgo (Lieu 
et al. 1996a), Coma (Lieu et al. 1996b), A1795 (Mittaz, Lieu and Lockman 1998) and 
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A2199 clusters (Lieu, Bonamente and Mittaz 1999). Recent attributions of these findings 
to spatial variations of the background (Bowyer, Berghoefer and Korpela 1999) were based 
upon an incorrect way of analyzing the EUVE data, and the decisive way of assessing 
the existence of CSE emission is that of an in situ background measurement by pointing 
time contiguously at small offsets from the cluster. This approach, applied to a recent 
EUVE re-observation of A2199 (Lieu et al. 1999b), clearly demonstrated that the offset 
pointing provides an accurate background template for the cluster observation, whereas an 
arbitrary blank field does not. With the aid of new EUVE observations, performed along 
with acquisition of in situ background, we were able to confirm the previous detection of 
CSE even in the most conservative case when we assumed no a priori understanding of the 
background distribution, so that a point-to-point subtraction of the background had to be 
adopted. We provide in this analysis two ways of utilizing the background data, one is to 
develop a model of the background and the other is direct subtraction of the offset field as 
in the case of A2199. 



2. Virgo 

A 46 ks EUVE observation of the Virgo cluster took place in February 1999, featuring 
a 34 ks pointing to the central galaxy M87 immediately followed by a 12 ks pointing at ~ 2 
degrees offset and with the same celestial orientation of the Lex/B filter. This observational 
strategy (see Lieu et al. 1999b for details) revealed the obvious detection of a halo of EUV 
emission around M87 of radius ~ 15' (Fig.l, top), confirming the original discovery (Lieu et 
al. 1996a). For both source and offset fields only the central, mildly vignetted parts of the 
DS Lex/B filter were considered ^ The in situ background measurement (Fig. 1, bottom) 

^Events with Y detector coordinate in the ~ 850-1200 range (Lieu et al. 1999c). This 
portion of the detector Y axis covers a linear size of ~ 26 arcmin; this implies that for radii 



-4- 



with the center of the annuh placed at the same detector position as that of the cluster's 
radial profile, indicates a distribution of surface brightness consistent with uniformity: 
the percentage deviation between the region occupied by the cluster and that beyond is 
statistically insignificant (the difference between the count rate of the 0'-17' and 17'-27' 
regions is 0.8±1.2 %). 

To account for any trends in the background radial profile which may be masked by 
statistical fluctuations of the measured data, we utilize the well estabhshed fact that the 
background consists of a spatially uniform particle component and a photon component 
which is vignetted (see section 4.2 and Fig. 7 of Sirk et al. 1997). Although the precise 
relative proportion varies with time, the former usually dominates, resulting in a radial 
proflle which is flat or slowly varying with no change in behavior at any particular radius, 
small scale ripples in the profile are negligible (Lieu et al. 1999c). Thus an empirical 
model of the background which involves a linear slope was satisfactorily fitted to the radial 
profile data at arcmin resolution. This model is then subtracted from the cluster data at 
same resolution with propagation of errors. Such a procedure was adopted, along with 
our constraints on the hot ICM parameters (see below) , to obtain the radial trend of the 
fractional EUV excess, the dashed diamonds in Fig. 2. 

For the more cautious reader, however, we also provided in Fig. 2 (dotted diamonds) 
the fractional EUV excess profile obtained from point-to-point subtraction of the in situ 
background and proper error propagation, when the subtraction was done after scaling away 
a small difference (~ 0.3%) in the absolute background level between the cluster and offset 
fields, by normahzing the 17' - 27' brightness of the two pointings (see Fig. 1). The small 
difference between the aforementioned background levels gives further indication that a 

>~ 13 arcmin the radial profiles are averaged only over those azimuthal angles that fall 
within the given Y axis detector limits. 
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contemporaneous, in situ offset pointing affords the most relevant background measurement 
(background levels in the DS detector have a dynamical range of ~ 300 %). Changes in the 
pattern of the background distribution is a higher order effect than changes in the absolute 
level of the background, as we demonstrated for our earher re-observation of A2199 (Fig. 
3a in Lieu et al. 1999b). It can be seen that even this approach, which yields the most 
conservative estimates of the soft excess, indicates positive signals and a rising radial trend. 
In fact, the two subtraction procedures gave statistically consistent results at all radii. 

In order to provide an accurate model of the hot ICM and its emissivity in the EUV 
band ^, an archival ROSAT PSPC observation of the Virgo cluster (RP number 800365) 
was analyzed. Here, effects of the gradient of Galactic line-of-sight HI column density 
(Nh) over the Virgo region, which are not entirely negligible, need to be evaluated. The 
Nff was shown by a recent 21 cm measurement (Lieu et al. 1996a) to radially increase 
from Nh = 1-8 x 10^° cm~^ at the cluster center to N^^ = 2.0 x lO^'^ cm~^ at a radius 
of ~ 17'. This was confirmed by IRAS 100 /im images (Wheelock et al. 1994; Snowden 
2000), and we consequently applied a higher Nh value to the outer regions of Virgo when 
modelling the spectral data. Given the radial HI variation (which continues its rising 
trend beyond 17'), and the known anti-correlation between HI and the 1/4 keV diffuse sky 
background (Snowden et al. 1998), one must assess how much the PSPC background was 
underestimated for the cluster region when it was determined, as in our case, from a ~ 
40'-50' annulus centered at M87. Of most concern is the 10'-15' annulus, where the 1/4 

^The PSPC data used in this paper were reduced according to the prescriptions of Snow- 
den et al. (1994) through the use of dedicated FTOOLS. Background fluxes are estimated 
from an outer annular region (near 45' from PSPC boresight) where the cluster contribution 
is considered negligible. Periods of high background (MV ratio > 170) were also rejected and 
PI channels 1-19 and 202-256, which suffer from calibration uncertainties, were not used. 
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keV sky background accounts for 12% of the detected flux in this band. A re-scahng of 
this background in accordance with the HI gradient (Snowden et al. 1998) over the annulus 
only leads to a negligible effect on the 1/4 keV excess (viz. a reduction by 1% from our 
reported precentage excess, see below). 

The emission from the hot ICM is at sufficiently low temperature (1~3 keV) to enable 
a spatially resolved temperature determination with PSPC data. PSPC R27 band X-ray 
spectra (PI channels 20-201, ~ 0.2-2.0 keV by photon energy; Snowden et al. 1994) were 
modelled, using the XSPEC software, with a photoelectrically absorbed single temperature 
MEKAL code (Mewe, Gronenschild and Van den Oord 1985; Mewe, Lemen and van den 
Oord 1986; Kaastra 1992), and the temperature and elemental abundances were fitted to 
the data. By means of the hot ICM parameters we assessed the existence of soft-excess 
in the EUV measurements. The results, as shown in Fig. 2, are positive. Further, they 
indicate the soft excess radial trend (SERT), or more precisely an increase in the fractional 
excess with radius. It should be mentioned that significant soft excess is also apparent in 
the C-band (synonimous with R2) as single temperature models do not fit well the PSPC 
data (Fig. 2, crosses; see also Lieu et al 1996a, Buote 2000). 



3. A1795 

A1795 was recently observed by EUVE for ~ 78 ks of exposure. The observational 
strategy utihzes the elongation of the DS Lexan/B filter (see Fig. 7 of Sirk et al. 1997). For 
about half of the observation, cluster emission was collected on one side of the detector, 
leaving the other side exposed only to background counts. For the remaining half of 
the time the role of the two sides was exchanged. Since our calculations based on the 
detected cluster EUV flux indicate that cluster's presence on the far side has no effect on 
the asymptotic background, this technique affords a time contiguous measurement of the 



- 7- 



cluster signals and the in situ background with equal exposures. The background radial 
profile was again found to be consistent with uniformity, as shown in Fig. 3, where the 
difference between the average count rate of the 2-11' and ll'-23' regions is 1.05 ± 0.97 
%. The CSE effect was detected out to a radial distance of ~ 10 arcmin from the center, 
confirming the previous EUVE observation (Mittaz, Lieu and Lockman 1998). 

Two archival ROSAT PSPC observations (RP numbers 800055 and 800105) were 
analyzed, to constrain the hot ICM parameters. The hne-of-sight HI is smoothly distributed 
at Nh = 1.05 xl0^° cm~^ in the direction of A1795 (Mittaz, Lieu and Lockman 1998), and 
the R27 band spectra of annular regions were fitted with a MEKAL model, the abundance 
being fixed at 0.31 solar (Fabian et al. 1994). The resulting EUV excess emission is shown 
in Fig. 4; like Virgo, it exhibits the SERT effect. We also found existence of CSE in the 
C-band, see Fig. 4 (crosses). 



4. Photoelectric absorption issues 

Throughout the analyses of this paper, we employed the absorption cross-sections of 
Morrison and McCammon (1982; MM82), whose He photoelectric cross-section is in good 
agreement with the compilation of Yan et al. (1998; Y98), the latter was employed in a 
recent investigation of the CSE (Arabadjis and Bregman 1999). In the C-band the MM82 
and Y98 He cross-sections differ by ~ 1%, whereas in the Lex/B passband (~ 65-190 eV) 
the Y98 cross-section is actually lower than in MM82 by ~8 %, i.e. the new cross-section 
renders our CSE fluxes even more prominent. Given that the cross section is not an issue, 
the proposition that the CSE detection of some clusters (including A1795 and Virgo) is in 
some way related to peculiarities in the Galactic absorption faces two problems. Firstly, the 
Nfj values needed to remove the CSE are grossly at odds with the ones measured by the 
dedicated observations of Mittaz, Lieu, & Lockman (1998). For example, we found that the 
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spectrum of the 7'-10' annulus of Virgo can satisfactorily be fitted by a one temperature 
plasma model (with no soft component) only by adopting an Nh as low as 1.34 ± 0.07 
xlO^° cm~^. Secondly the SERT effect (e.g., Fig. 2 and Fig. 4) implies a curious spatial 
correlation between our Galactic interstellar medium and the centers of the two clusters. 
We therefore conclude that the CSE is a genuinely celestial emission component. 

5. Interpretation of the soft component and conclusions 

We now discuss the physical significance of the soft component. Currently there are 
two contesting models: the non-thermal scenario, which invokes an Inverse Compton (IC) 
effect as cause for the emission (Sarazin and Lieu, 1998; Hwang 1997; Ensslin and Biermann 
1998), and the original thermal model, which explains the CSE as thermal emission from a 
'warm' gas (~ 10^ K; Lieu et al. 1996a). 

A recent study of the Coma cluster (Lieu at al. 1999a) showed that the CSE emission 
could be explained as an IC effect, albeit invoking a large population of energetic electrons. 
When the present Virgo soft excess was fitted with a power-law, we also accomodated the 
possibihty of 'aging' due to loss of the relativistic electrons to the ICM. ^ Even so, the fit 
was far from acceptable. A simple power-law model for the soft component requires an 
extremely large value of the differential photon number index (Table 1), inconsistent with 
the value of ~ 1.75-2 expected from a population of electrons accelerated at astrophysical 

^A spectral code which introduces radiative (synchrotron and IC) and non-radiative 
(Coulomb) losses on a power-law model was integrated with the XSPEC spectral fitting 
package. As function of 'aging' time {iage) and environmental conditions (such as density of 
the hot gas and magnetic field strength), electrons at all energies are progressively removed, 
modifying the resulting emissivity. 



-9- 



shocks (see Lieu et al. 1999a and references therein). When aging effects are included, 
a satisfactory fit can be found only for some of the regions, and the model was rejected 
for regions of strongest signal (0'-7' from M87; see Table 1). The CSE can be explained, 
however, as thermal emission from a thin plasma at lower temperature (~ 0.1 keV, see 
Table 1), following the original interpretation of Lieu et al. (1996a). 

Fitting of the CSE of A1795 with a non-thermal model produced a similar outcome. At 
large radii, the EUV excess is accompanied by a weaker C-band excess. This necessitates a 
sharp cut-off, which can be understood as due to aging, and a highly evolved IC spectrum 
is a possible model. Table 2. However, the energetic requirements are severe, see Table 3, 
because the electron pressure alone (in the absence of relativistic ions, which accentuates 
the difficulty ~ 10 - 100 times) often exceeds that of the hot ICM. The 6-10 arcmin region 
requires relativistic electrons which, at time of injection, had a pressure ~ four times higher 
than that of the hot gas, leading to major confinement problems. The thermal scenario 
again appears more appropriate: the data of Table 2 indicate that a ~ 0.05-0.1 keV gas 
can satisfactorily explain the CSE. 

What are, then, the implications of the thermal interpretation on cluster masses? In 
absence of any knowledge on the clumping of the warm phase, a 100 % filling factor implies 
a warm gas mass comparable with that of the hot gas (Table 4). This is in good overall 
agreement with recent large-scale hydrodynamical simulations (Cen and Ostriker 1999) 
which indicate that at the present epoch the baryons are to be found, in similar amounts, 
at hot (~ 10® K) and warm (~ 10^ — 10^ K) temperatures, and only marginally at lower 
temperatures (see also Bonamente, Lieu and Mittaz 2000 for the detection of cold clouds 
with commensurately smaller mass budgets). 
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Figure captions 

Figure 1: Radial profile of the cluster pointing (top) and offset pointing (bottom) of 
the Virgo observation. For both plots the radius is measured from the detector position of 
M87. The innermost 1' bin of the cluster, dominated by the emission from M87, is omitted. 
The dotted line represents average brightness of the 17'-27' region. On same ordinate of 
the cluster pointing (top) is overplotted the radial profile of the offset pointing with open 
circles, after a ~ 0.3 % rescahng (the difference between the 17'-27' region brightness of the 
two profiles) is applied. A radial profile of offset pointing (bottom) in 1 arcmin resolution 
was successfully fitted to an empirical linear slope, and its best-fit parameters used for the 
evaluation of the fractional EUV excess (dashed diamonds in Fig. 2; see also text). 

Figure 2: Diamonds are radial profile of the EUV fractional excess r] of Virgo, defined 
as Tj — (p — q)/q, where p is the observed EUV band flux after background subtraction, 
and q is the expected flux from the hot ICM. Vertical semi-diameters are here and after 
1-cr errors. Dashed diamonds are obtained by subtraction of model background, dotted 
diamonds by point to point background subtraction. Crosses are radial proflle of the C-band 
fractional excess rj of Virgo. 

Figure 3: Radial proflle of the co-added cluster pointings (top) and offset pointings 
(bottom) of the A1795 observation. The innermost 2' bins are omitted, as there was 
substantial contamination from a bright unidentified transient source near the cluster 
center. On either plot the dotted line represents average brightness of the 11 '-23' region. 
On same ordinate of the cluster pointing (top) is overplotted the radial profile of the 
offset pointing with open circles, after a ~ 0.7 % rescahng (the difference between the 
11 '-23' region brightness of the two profiles) is applied. A radial profile of offset pointing 
(bottom) in 1' resolution was successfully fitted to an empirical linear slope, and its best- fit 
parameters used for the evaluation of the fractional EUV excess (dashed diamonds in Fig. 
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4). 

Figure 4: Diamonds are radial profile of the EUV fractional excess r] for A 1795. Owing 
to contamination from a transient source (see Figure 3), data for the 0-2 arcmin bins are 
from a previous observation of the cluster (Mittaz, Lieu and Lockman 1998). Note also 
that this previous paper showed a similar plot with a stronger rising trend, because it used 
absorption code of Balucinska- Church and McCammon (1992), whereas the present work 
uses the original Morrison & McCammon (1982) code, which is now believed to involve a 
more reliable He cross-section (see text). Dashed diamonds are obtained by subtraction of 
model background, dotted diamonds by point to point background subtraction (see section 
2). Crosses are radial profile of C-band fractional excess of A1795. 
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Table 1: Modelling of Virgo spectra with a two-temperature MEKAL code (left), MEKAL 
+ power-law model (center) and MEKAL + 'aged' power-law model (right). Errors are 90 
% confidence (x^ -|- 2.701 criterion). For the latter, the IC magnetic field was fixed at 1//G, 
the density of the ICM calculated from the best-fit /5-model of Nulsen and Bohringer (1995) 
and the electron differential number index 7 fixed at the value of 2.5 (wich corresponds to a 
similar photon index a of 1.75). The parameters of the hot phase were fixed at those best 
values obtained by fitting PI channels 42 - 201 with a single temperature model of floating 
abundance. Nwarm is here and after in units of 10~^^/47rD^n^V, where D (cm) is the distance 
to the source, n (cm~^) is the gas density and V {cvci?) is the volume of the emitting region. 
For MEKAL-FMEKAL and MEKAL-^PO model fits, reduced is always between 1.1 and 
1.3 for 180-181 d.o.f.. 
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Table 2: Modelling of A1795 spectra with a two-temperature MEKAL code (left) and with a 
MEKAL + 'aged' power-law code (right); errors are 90% confidence (x^ + 2.701 criterion). 
Parameters of the hot phase were fixed at those best values obtained by fitting PI channels 
42 - 201 with a single temperature model of abundance 0.31 solar (Fabian et al 1994). For 
the non-thermal model, the value of the IC magnetic field was fixed at 1/iG; the density of 
the ICM, a parameter which affects the aging of electrons, was calculated from the best-fit 
/3-model of Brief and Henry (1996). Electron differential index 7 was fixed at the value of 
2.5 (see Table 1). 

MEKAL+MEKAL MEKAL+'aged' PO 

Region T^arm. Nu,arm red. tage Red. 

(arcmin) (keV) xlO^ (d.o.f.) (Gyrs) (d.o.f.) 

1- 2 0.125±g:gi| 0.045±g:gg| 1.0(363) 1.7 ig;}! 1.0(363) 

2- 3 0.11±g:gf 0.0285±g:gg? 1.16(363) 2.25 ±111 1.15(363) 

3- 6 0.049±[]:[]}^ 0.12±g;^^ 0.97(363) 2.8 0.98(363) 
6-10 0.025±[]:g}| llOib^go 0.88(363) 3.7 ±lf^ 0.88(363) 
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Table 3: Luminosity and pressure estimates for A 179 5; L](2 intrinsic luminosity of 

non-thermal best-fit model between photon energies 65-250 eV (7mm=279 and 7maa;=548) 
in units of 10^^ (erg s~^) and P^'se jg ^j^g pressure of relativistic electrons in the [^minilmax] 
Lorentz factor interval. and P'^{t = 0) are the total electron pressure at the present epoch 
and at the acceleration epoch, respectively, assuming cosmological parameters Hq = 75 km 
s~^ Mpc~^ and = 0. If Hq — 50 km s~^ Mpc~^ is adopted, all present pressure estimates 
for the non-thermal component will increase ~ two-fold. Pressure for the ICM (ps"^) was 
calculated according to the parameters of Brief and Henry (1996). 
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Table 4: Warm gas mass estimates and comparison with the hot ICM for the Virgo and 
A1795 cluster (see Tables 1 and 2). Warm gas densities (n^ar-m) are estimated for 100% 
filling factors. 
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